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Abstract 



We discuss perturbative and non-perturbative strong interaction effects in the pair 

. ^ production of stop squarks (ti) at e+e~ colliders. Events with an additional hard gluon 

Q.' allow to detect or exclude titj production even in scenarios with very small mass splitting 

between ti and an invisible lightest supersymmetric particle (LSP). Such events can also 

Qj \ help to establish that ti transforms as a triplet under SU{3)c- We also carefully study 

i-G ' non-perturbative ti fragmentation, which is currently not well understood: not only is 

the ti fragmentation function not known very well, but also there are ambiguities in 

k><( I the algorithm employed to model fragmentation. We present numerical results both for 

^ ' CERN LEP-183 and for a proposed future e+e~ collider operating at center-of-mass 

energy -y/i = 500 GeV. 



1) Introduction 



The lighter scalar top (stop) mass eigenstate ti is likely to be the lightest squark ^. In 
models where supersymmetry breaking is communicated to the superpartners of the known 
SM particles at some energy well above the weak scale, renormalization group effects reduce the 
soft mass of the stops compared to those of other squarks. In addition, after the spontaneous 
breakdown of the electroweak SU{2) x f/(l)y symmetry, SU{2) doublet and singlet stops mix, 
leading to a further reduction of the mass m^^ of the lighter eigenstate. In fact, ti might even 
be the lightest charged sparticle. It is therefore actively being searched for at both pp and 
e"*"e~ [^, ^ colliders. In this article we focus on iii\ production at existing and future e"*"e~ 
colliders. 

Unlike the other candidates for the lightest charged superparticle, sleptons and the lighter 
chargino mass eigenstate xf ; ii has strong interactions. Stop pairs can thus be produced 
in association with a hard gluon. Events of this type can help to overcome Standard Model 
backgrounds to stop production in scenarios where ti decays into a lightest neutralino x5 ^^^ ^ 
quark, if x\ is invisible (stable or long-lived) and mi^ — rriy. is small. In such cases the emission 
of a hard gluon can increase the total visible energy as well as the total missing transverse 
momentum, helping to discriminate between titl events and backgrounds from two-photon 
reactions. Once ti has been discovered, one will also want to study its properties. Events with 
an additional hard gluon will then be crucial in establishing the color charge of ti. 

If m^^ < mh + m-± , rrit + m^, ti is expected to be quite long-lived . It will then fragment 
into a "stop meson", which is a spin-1/2 particle, plus a fragmentation jet prior to its decay. 
This is a non-perturbative process, and hence not well understood. In fact, the very concept 
of fragmentation is somewhat ill-defined for massive particles: whereas one massless parton 
can fragment into two collinear massless partons while conserving energy and momentum, 
a massive particle cannot fragment into another particle with equal (or larger) mass and a 
second, massless particle with nonvanishing energy. This does not pose much of a problem in 
the calculation of single-particle inclusive cross sections, through which fragmentation func- 
tions are usually defined, but it leads to ambiguities in the implementation of fragmentation 
in event generator programs. This is of special importance in scenarios with small ti-x? mass 
splitting, where the fragmentation jets may contribute the bulk of the visible energy. 

Since non-perturbative ti fragmentation occurs even if no hard gluon is emitted, we begin 
by discussing this process in the simple case of e+e~ -^ tii^. We introduce three different 
fragmentation schemes, which differ in the precise definition of the fragmentation (scaling) 
variable x, and/or in the algorithm used to ensure overall energy-momentum conservation. 
These three algorithms lead to quite different stop energy distributions after fragmentation 
when used naively with the same fragmentation function. On the other hand, these differences 
become quite small if the parameter of the fragmentation function is adjusted such that the 
mean stop "meson" energy is the same in all three scenarios. 

In Sec. 3 we discuss events where one hard gluon is emitted. We give an explicit expression 
for the differential cross section in a form that is easily implementable in an event generator. 
We then demonstrate, using a parton-level simulation, that the possibility of emitting a hard 
gluon allows to discover or rule out ti even for very small mass splitting to the LSP. To this 
end we employ a set of cuts modeled on the OPAL analysis 0] . These cuts also greatly reduce 
the sensitivity to the details of ti fragmentation. We also present a second analysis for a 
future e+e~ collider operating at center-of-mass energy -^i = 500 GeV, as an example of the 
type of analysis one might perform to determine the color charge of ti. In this case the cross 



section after cuts does depend on the details of fragmentation, but the resulting ambiguity in 
the prediction is much smaller than the perturbatively calculable contribution. Finally, Sec. 4 
is devoted to a brief summary and conclusions. 



2) Fragmentation 



We begin our analysis with a discussion of ti fragmentation. We will model it using the 
Peterson et al. fragmentation function 0: 
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where the single free parameter e^ is expected to scale as the inverse of the square of the stop 
mass 



q=— , (2) 

tl 



with e ~ 0.1 to 0.5 GeV^. Di has to satisfy the "stop number conservation" condition 
/^ . dx xDi{x) = 1, where Xmm is the kinematic minimum of x. This leads to 

.,(...„.,, . _^ _ 1 ,„, (i^ , 1 _ _J_ ,0iV^>^ (3) 

Fragmentation functions are usually defined through differential single-particle inclusive 
cross sections. In the simple case of leading order stop pair production, e~^e~ -^ iiil, the stop 
energy prior to fragmentation is fixedQ to -\/i/2, so that 

Unfortunately the choice of Q, and consequently of the associated fragmentation variable x, 
is somewhat ambiguous for massive particles. It can be the energy E of the stop "meson" tM 
after fragmentation, in which case 

x = XE = ^ = ^ ; (5) 

or one could choose the absolute value of the 3-momentum k, 

x = Xk = -^ — = —r-F- ' (6) 
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where (3 = Jl — 4m? /s. The second identities in eqs. (|) and (||) hold only for the simple 
case of 2 -^ 2 kinematics. Occasionally the fragmentation variable is also defined through the 
"light cone variable" {E + \k\)/2. This will give results intermediate between the other two 
choices, so we focus on definitions (^) and @ here. 

*We ignore initial state radiation in our analysis. It can reduce the total cross section by ~ 20% R], but 
will not change most kinematical distributions very much, which are the main emphasis of this work. 



A second ambiguity is related to the fact that a massive (s) quark cannot fragment into a 
hadron of equal or greater mass and one (or more) massless parton(s) without violating energy 
and/or momentum conservation, if all particles are on-shell. It is therefore not obvious how 
to introduce fragmentation functions, defined through inclusive cross sections as in eq. (^, 
into event generator programs that attempt to describe all (relevant) properties of events. 
We have implemented two different classes of algorithms that ensure energy and momentum 
conservation at the level of the global event. For simplicity we always model the "fragmentation 
jet" as a single massless parton collinear with the parent stop. 

The first class of algorithms is based on rescaling the 3-momenta of all particles in the final 
state by some overall factor, after the stops have been fragmented. The value of this factor is 
calculated for each event through iteration, by requiring that all final particles are on-shell and 
that their energies in the e~^e~ center-of-mass (cms) frame add up to a/s. These algorithms 
require that the 3-momenta of all particles after fragmentation, but before rescaling, add up 
to zero (in the e"*"e~ cms frame); overall 3-momentum conservation is then ensured throughout 
the rescaling process. This means that the fragmentation step itself always has to fix the 3- 
momentum of the fragmentation jet, even if the fragmentation variable is defined through the 
stop energy as in eq. (^). In other words, this fragmentation prescription "locally" conserves 
3-momentum, but violates energy. Overall ("global") energy conservation is then restored 
through the rescaling procedure, which not only changes (reduces) the 3-momenta of all 
particles in the final state but also gives rise to a flow of energy between different parts of 
the event. This class of algorithms is similar to that employed for the fragmentation of heavy 
quarks in the ISA JET event simulation package ||^. It can be combined with either definition 
(^ or (^ of the fragmentation variable; we will call the two resulting algorithms Ie and I^, 
respectively. 

Note that the sum of the square of the 4-momenta of a stop "meson" and its associated 
fragmentation jet after rescaling exceeds m? in the prescriptions Ie and Ik- In other words, one 
implicitly assumes that the stop squarks were produced slightly off-shell when implementing 
fragmentation.^ Our second algorithm makes this explicit, by directly generating a distribution 
in the squared (off-shell) mass m^ of the produced stop squarks: 
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where x is distributed according to the fragmentation function -Dj(x) of eq. (|I|), and E^_^ is the 
energy of on-shell stops, as generated during the calculation of the stop production matrix 
element. The form (|^) has been obtained by writing the stop 4-momenta prior to and after 
fragmentation as {E^^,k) and {xEf_^,k'), respectively, with k' being parallel to k. One then 
requires that the difference of these two 4-momenta is light-like (for a massless fragmentation 
jet), and that the stop meson has mass m^^. This allows to express m^ as a function of 
Ef_^ and k. In the last step one replaces these quantities by the energy and 3-momentum of 
on-shell ti; this can be regarded as a one-step iteration. Had we instead defined the stop 
"meson" 4-momentum as {E',xk), i.e. tried to implement definition (P) in this scheme via a 
one-step iteration, the second term on the r.h.s. of eq. (0) would have to be multiplied by an 

TThe matrix elements for stop production are nevertheless computed assuming ii and i^ to be on-shell, 
so as to reproduce the correct inclusive titl production cross section, in accordance with QCD factorization 
theorems. The same apparent contradiction is built into all QCD fragmentation or parton shower programs 
we know of. 



additional factor Jl — m~ /E? . However, we will stick to the ansatz (|^) for scheme II here, 

since the difference between the definitions (^ and (^) of the scaling variable will be explored 
by comparing the results of schemes I^; and I^. 

Once the masses of ti and il have been generated according to eq. (|^, the energies and 
3-momenta of all particles in the final state have to be recomputed, using new limits of phase 
space as appropriate for the increased stop masses. In so doing, care has to be taken to 
generate kinematic quantities that "correspond" to those used (prior to fragmentation) to 
compute the squared matrix element for stop production, which enters the event weight. This 
can be accomplished as follows. Assume the phase space is parameterized by a set of kinematic 
quantities g^, i = 1, . . . ,N. The kinematic limits of g^ are determined by ^/s, by the masses of ti 
and tl, and (for i ^ 1) by the values of gi, . . . , qi-i. An event can then be characterized by a set 
of dimensionless numbers rj G [0, 1], so that qi = gi,min + ^j(gj,max — gj,mm); the r^ can be viewed 
as the random numbers used by a Monte Carlo integration of phase space. The correspondence 
between "fragmented" and "unfragmented" events can then be ensured by using the same set 
{rj} when regenerating all 4-momenta after fragmentation, only modifying the kinematical 
limits to accommodate increased stop squark masses for the "fragmented" events. Once the 
4-momenta of the off-shell squarks have been generated in this way, fragmentation itself can 
be treated as the (coUinear) decay of a particle with mass m into a second particle with mass 
m^^ and a massless parton. 

This second algorithm has both advantages and disadvantages. On the positive side, there 
is no need for iterative rescaling of 3-momenta after fragmentation; however, as described 
in the previous paragraph, all 4-momenta do need to be recomputed, allowing for different 
(virtual) ti and il masses. Another potential advantage is that, although the ansatz (|^) 
has been designed to emulate a conventional fragmentation function using definition (^ of 
the scaling variable, it could easily be converted into a fully Lorentz invariant fragmentation 
scheme, since the mass m^ is a Lorentz invariant. This would be more difficult to accomplish in 
both variants of scheme I, which are based on non Lorentz covariant quantities (energies and 3- 
momenta). On the other hand, at least as implemented using eq. (|^), scheme II does not have 
a smooth limit as m^^ — > 0. Eq. (|^ then always gives m.^ = 0, independent of x, leaving the 4- 
momenta of the fragmentation products undefined. In contrast, in the massless limit schemes 
I^; and Ifc both approach the same scheme (corresponding to independent fragmentation, again 
as in IS A JET ^), where no rescaling of 3-momenta is required. 

Rather than working with variable (virtual) stop squark masses and fixed stop "meson" 
mass, one might also consider keeping the stop squark masses fixed at m^^ , but varying the stop 
meson masses. This would also allow to implement energy-momentum conservation without 
any rescaling. However, there are several problems with such a scheme. Energy-momentum 
conservation would force the "mesons" to be lighter than the squarks, whereas fiavored D and 
B mesons are heavier than c and b quarks. Also, a continuous stop meson energy distribution 
would correspond to a continuous stop meson mass spectrum as well. For these reasons we 
will not explore this alternative any further. 

In Figs. la,b we compare results for our three fragmentation schemes. We have chosen a 
parameter set of current interest, ^/s = 183 GeV and mj^ = 80 GeV. The solid, long dashed 
and short dashed curves refer to schemes I^, 1e and II, respectively. We have taken ti to be 
a pure SU{2) doublet {6i = in the notation of refs. [|7|, ^), but this only affects the overall 
normalization of the curves. We have also included "soft and virtual" QCD corrections 0, 
defined through an infrared cut-off for emitted gluons of 1 GeV. This ensures that we are only 
studying events with 2 — > 2 kinematics; events with additional hard gluon will be discussed 



in Sec. 3. 

In Fig. la we compare the stop "meson" spectra produced by our three schemes. We 
have taken the same value of e in all three cases, but the resulting spectra clearly differ 
substantially. The difference between schemes Ie and Ifc is fairly easy to understand. Since 
the produced stop squarks are far from being ultra-relativistic (/? ~ 0.5), reducing their 3- 
momenta by some factor x corresponds to a much smaller relative reduction of their total 
energies, which receive a large (constant) contribution from their mass. As a result, for fixed 
e scheme I^ leads to a much harder stop "meson" spectrum than scheme Ie does. Except 
for the steep drop towards the lower kinematical endpoint of the distribution, scheme II falls 
in between the two variants of scheme I here. The large differences displayed in this figure 
clearly show that care must be taken when comparing the values of the free parameters of 
heavy (s) quark fragmentation functions obtained from different event simulation programs, 
with possibly different implementations of the fragmentation process. 

On the other hand. Fig. lb shows that at least for the simple case of exclusive titl pro- 
duction (2 — i> 2 kinematics), the differences between fragmentation schemes can largely be 
absorbed into a change of the value of the parameter e of eq. (0). Here we have chosen e 
in each scheme such that the average stop meson energy coincides with the prediction of the 
original Peterson et al. function (|l|); using e = 0.5 GeV^ and the definition (|]) of the scaling 
variable, it predicts (-E^^^) = 89.1 GeV. Almost no change of e is required in scheme II. On 
the other hand, in scheme I^;, e has to be reduced to 0.1 GeV^, while in scheme I^ it has to 
be increased to 2.0 GeV^, in order to reproduce the same {Ei^^). These two values of e differ 
by a factor of 20! 

After this adaptation of e, the three schemes give quite similar results for the stop "meson" 
spectrum, except near its lower kinematical boundary. However, recall that these schemes 
ensure energy-momentum conservation only globally. This means that the fragmentation 
process will change the 4-momenta of partons that seemingly have nothing to do with it. This 
is actually not as unreasonable as it may sound, since color neutralization in any case connects 
different partons during hadronization. Furthermore, notice that none of the three schemes 
reproduces the original Peterson function, shown as the dotted curve, near the upper endpoint 
of the spectrum. The reason is that, due to the requirements of overall energy-momentum 
conservation, our schemes allow Ei^^ ~ -\/^/2 only if both fragmentation variables are close to 
1, which is obviously much less likely than having only one of them close to unity. We will see 
in the next Section that even after adapting the value of e so as to produce the same {E^^^), 
the three schemes lead to somewhat more pronounced differences in events with additional 
hard gluon. 

3) Hard gluon emission 

In this Section we discuss the production of a tit^ pair in association with a hard gluon, 
e^e~ -^ tit\g. The gluon can be emitted from either of the two stop lines; in addition, 
there are diagrams with a 'ytitlg or Ztitlg vertex. This process has been treated previously in 
refs. [^ [To], |l^ , using parameterizations of phase space that simplify analytical calculations. In 



contrast, our treatment is geared towards easy implementation in an event generator program. 
The spin-averaged squared matrix element is 
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Figure 1: The stop "meson" spectrum as predicted by tlie three fragmentation schemes de- 
scribed in the text. All three schemes employ the Peterson et al. fragmentation function. In 
a) the same value of the parameter e has been used in all cases, while in b) the value of e has 
been adjusted so as to reproduce the prediction (E^^^) ~ 89 GeV of the original Peterson et 
al. fragmentation function with e = 0.5 GeV^, using the definition (|^) for the scaling variable; 
this gives e = 0.1,2.0 and 0.5 GeV^ for schemes Ig, 1^ and II, respectively. The shape of the 
Peterson et al. fragmentation function is shown by the dotted curve. 
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where Qs is the strong couphng constant, which we take at scale /i = ^/s, and e is the QED 
couphng. The coefficients A and B describe the strength of the effective vector and axial 
vector couplings, being given by 



A = qi- 
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describes the Ztitl coupling [0 with 9^ being the stop mixing angle (the angle between ti and 
^l)', qi (= +2/3) and Is- (= +1/2) are the ti charge and t^ weak isospin, respectively. 

The cross section can now be expressed in terms of the reduced squared matrix element 



M 



introduced in eq. (||): 
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Here, Eg and Ei are the energies of the gluon and of ti in the e^e~ cms frame, 9g is the angle 
between the gluon and the direction of the positron beam, and 0j determines the direction of 
the ti 3-momentum, as follows. We choose the frame where the gluon's 3-momentum is given 
by 

g'=E,(0,sine„cos^J , (12) 

i.e. the integration over the polar angle of the gluon has been performed explicitly in eq. (|TTp, 
yielding a factor of 27r. The ti 3-momentum is then given by 
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2Eg^Ef~ml 
is the angle between ti and the gluon, which is fixed for given Eg and E^:. Note that we have 



allowed for different ti and tl masses mi and m2 in eq. (14), as required for the implementation 
of our fragmentation scheme II discussed in the previous section. However, for the purpose 



of computing event weights rrii = m2 = mi^ should be used. The t\ 3-monientuni can now 
trivially be computed from eqs. (p!2D-(p^) as ^2 = —ki — q. 

Finally, the limits of phase space integration, again allowing for different ti and t*^ masses, 
are: 
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Note that the squared matrix element (|^) diverges linearly as Eg —* 0. This requires the 
introduction of an infrared cut-off ii^™™ in eq. (|15a|) . The same cut-off also appears (logarith- 
mically) in the expression for the "soft + virtual" QCD corrections. Physical quantities do 
not depend on -E™™, if it is chosen sufficiently small so that events containing a gluon with 
energy Eg = E"^^^ still effectively possess 2^2 kinematics. We have checked numerically 
that eqs. (3)-(|T6|) reproduce the "hard" real QCD corrections to the total iiil production cross 
section as given in ref. (and contained implicitly in |^). Of course, stop fragmentation, as 
implemented using either of our schemes, does not change the total cross section before cuts. 

In the subsequent analysis we have taken into account the ti decay, which we have assumed 
to be always into a massless quarkQplus an invisible LSP Xi- A first illustration of the impact of 
both perturbative and non-perturbative strong interaction effects on tit* production is shown 
in Fig. 2, where we plot the total visible energy for y/s = 183 GeV, m^^ = 70 GeV, and m^ = 67 
GeV. We have chosen a scenario with very small ti-LSP mass splitting, since in this case strong 
interaction effects are most important. This is illustrated by the difference between the dotted 
curve, where all such effects have been switched off, and the other curves. In the absence of 
strong interaction effects, and before any cuts have been applied, the visible energy spectrum 
has a triangular shape, centered at E* ■ (v^/mjj ^ 7.7 GeV, where E* = (m? — m?)/(2mfj 
is the quark energy in the ti rest frame. We see that strong interaction effects are expected 
to nearly double the average visible energy for our choice of sparticle masses. 

However, there are significant uncertainties in this prediction, due to our lack of under- 
standing of non-perturbative ti fragmentation. Different fragmentation schemes using the 
same fragmentation function can give quite different results, as can be seen by comparing the 
long and short dashed curves, which are valid for schemes Ie and II, respectively, with e = 0.5 
GeV^. Such a difference is not surprising, given the results of Fig. 1. However, even after 
adjusting the value of e used in scheme Ie so as to reproduce the average stop "meson" energy 
predicted by scheme II, which is very close to the value predicted directly from the Peterson 



*In models with exact universality of the squark masses at the scale Mx of Grand Unification one expects 
ii to dominantly decay into a charm quark ||^; in this case the visible energy could be further degraded 
by semileptonic decays of the charm (anti)quark. However, there is no compelling reason to assume exact 
(rather than only approximate) squark degeneracy at Mx- Moreover, ii decays are independent of the strong 
interaction effects on tit* pair production, which are the main focus of our analysis. 
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Figure 2: The total visible energy in e"'"e~ -^ titl{g) events at LEP-183. The dotted curve 
has been obtained by switching off both perturbative and non-perturbative strong interaction 
effects, while the dot-dashed curve includes non-perturbative ti fragmentation but no hard 
gluon radiation (defined with an infrared cut-off i?™™ = 1 GeV). The other curves include 
both perturbative and non-perturbative effects. The labels refer to the fragmentation scheme 
used, and give the value of the fragmentation parameter e of eq. (0) in units of GeV^. 



et al. fragmentation function, some discrepancies remain. Indeed, the difference between the 
predictions of scheme \e with adjusted e parameter (solid curve) and scheme 11 is comparable 
in size to the effect of hard gluon radiation, which can be seen by comparing the short dashed 
and dot-dashed curves; the former shows the full cross section, while the latter is for the 2 — > 2 
contribution only. 

One might be tempted to conclude from Fig. 2 that non-perturbative stop fragmentation by 
itself can release a sufficient amount of visible energy to give a viable signal for tit\ production 
even for small ti-LSP mass splittings. However, one should keep in mind that in the absence 
of hard gluon radiation, the fragmentation products will always be almost back-to-back in 
the transverse plane, and the ti decay products are not energetic enough to change this 
picture. Unfortunately there are large backgrounds to such configurations, even if one requires 
a substantial transverse momentum imbalance. The main source of backgrounds are two 
photon processes. They can lead to large (apparent) missing pt either through semileptonic 
decays of charm quarks, or through mismeasurement since the energies of the jets are quite low, 
so that calorimetry becomes rather imprecise. Recent experimental searches |^, Q therefore 
discard such back-to-back configurations. This will enhance the importance of signal events 
containing an additional hard gluon. 

This is illustrated in Fig. 3, where we have attempted to incorporate the most important 
cuts applied by the OPAL collaboration p[ into our parton level event generator. We show the 
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Figure 3: The total inclusive t\t\ cross section at LEP-183 after the cuts ( p!7D have been 
applied. The solid curves show the full result, while the dashed curves have been computed 
by switching off hard gluon radiation. For mi^ = 50 and 60 GeV, the infrared cut-off has 
been chosen such that the total cross section before cuts of the 2^2 contribution reproduces 
the leading order, C(a°), prediction; this gives i?™™ = 5.4 (1.8) GeV for m^^ = 50 (60) GeV. 
For ?Tij^ = 70 GeV this would give a very small value of E^ 
instead. Recall that the total result (solid curves) does not depend on Eg 
for fragmentation scheme 11, except the dotted one, which is for scheme 1^. 



™^'^, so we have chosen 1.0 GeV 
™"^ All curves are 



total inclusive titl cross section at i/i = 183 GeV after the following cuts have been applied: 

M^i, > 4.0 GeV ; (17a) 

PT,miss > 4.5 GeV ; (17b) 

(17c) 
>20° . (17d) 



cos6'miss| < 0.9 ; 



^acop 



In order to model the geometric acceptance of a typical LEP detector, we have discarded all 
partons with azimuthal angle | cos^| > 0.97 before computing the total visible mass M^is and 
missing transverse momentum PT.miss- The cut ( |17c|) is on the angle of the missing momentum 



vector with the beam pipe. OPAL also requires that the kx clustering algorithm [|12[ with 
fixed jet resolution parameter j/cut produce exactly two jets, and applies cuts on the invariant 
masses of these jets. However, it is difficult to model such cuts accurately with a parton level 
event generator. Furthermore, insisting on exactly two jets with given j/cut would tend to 
remove events containing a hard gluon in addition to the titl pair. Therefore, the full set of 
cuts applied by the OPAL collaboration is not efficient for scenarios with very small ti — Xi 
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mass splitting. We have instead used the /c^ clustering algorithm to force the event into a 
two-jet structure, except for events with less than two partons in the acceptance region of the 
detector, which we discarded - that is, we kept merging jets until exactly two jets remain. We 
then applied the cut ( |17d| ) on the angle between these two jets in the transverse plane. 



The results of Fig. 3 indicate that the emission of a hard gluon should allow to detect or 
exclude titl pair production even for arbitrarily small ti — x? mass splitting. The published 
OPAL limit 0] on m^^ corresponds to a total cross section of 0.13 pb, after the cuts (p!7|) 
have been applied. This indicates that existing data might allow to exclude, or find evidence 
for, a ti nearly degenerate with the LSP up to m^^ ^ 60 GeV for mixing angle 6^ = 0. Of 
course, a complete experimental analysis would be needed to really establish such a bound. 
We emphasize that an accurate description of hard gluon radiation would be crucial for such 
an analysis. Leading-log parton shower models, as used in refs. 0, ^, cannot be expected to 
model this process accurately. The importance of this 2 — > 3 contribution is illustrated by the 
dashed curves, where this contribution has been turned off; the cross section after cuts then 
goes to zero in the limit of vanishing ti — x? niass difference. 

Most results shown in Fig. 3 have been computed using scheme II with e = 0.5 GeV^. 
However, the cuts ( |17|) greatly reduce the sensitivity to the details of the fragmentation process. 
This is illustrated by the dotted curve, which shows the inclusive tit* cross section after cuts for 
m^-^ = 70 GeV, as predicted using scheme I^ with the same value of e. We saw in Fig. 1 that this 
corresponds to a significantly harder stop "meson" spectrum, i.e. much softer fragmentation 
jets. Nevertheless the predicted cross section is almost the same as that for scheme II. This 
is due to a somewhat fortuitous cancellation: models with harder fragmentation jets predict 
a significantly larger cross sections after the cuts ([T7|a,b) have been applied; however, the 
resulting event then also tends to be more back-to-back, so that a relatively larger number of 
events is removed by the cut ( |1 7d|) . This weak sensitivity to the details of the fragmentation 



process has also been noticed in the experimental searches P, ^, where it was studied by 
varying the parameter e of the Peterson et al. fragmentation function. This result is very 
welcome, since it almost removes one (potentially large) source of uncertainty of the measured 
(limit on the) total tit* cross section. 

As already mentioned in the Introduction, events with an additional hard gluon might also 
help to determine the color charge of ti.[|] Given the present experimental limits, a definitive 
analysis of this kind may only be possible at a high energy, high luminosity linear e~^e~ collider 
(LC)I^ Various designs of such colliders are now being studied, with quoted luminosities ranging 
up to 500 fb-i/yr g|. 

Assuming that ti always undergoes 2-body decay into a light quark and an invisible LSP 
Xi, the most direct way to determine its color charge is probably through the study of events 
with 3 or more jets. Notwithstanding, non-perturbative ti fragmentation can also produce jets. 
Fortunately these fragmentation jets are expected to be quite soft; perturbative contributions 
should therefore dominate for sufficiently large jet energies. 

This is borne out by Fig. 4, which shows results of an analysis for ^/s = 500 GeV, tti^^ = 150 
GeV, and rriy. = 100 GeV. Allowing for the (probably) larger beam holes of future LC detectors 

^The "3-jet" contribution to the total inclusive cross section for the production of a pair of colored scalars 
has been calculated in |1^; however, that analysis did not treat fragmentation and decay of the scalars, nor 
were any cuts applied to reject the background. 

^Oi course, the discovery of ii at a hadron collider would immediately prove that it carries a nonvanishing 
color charge. However, in order to confirm that it belongs to an SU{3)c triplet one would have to measure 
both m^^ and the iiij production cross section, which may not be easy. 
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as compared to present LEP detectors, we have discarded all partons with | cos6\ > 0.9. We 
then used the fc^ clustering algorithm [Q to force the event into a 3-jet structure. Following 
ref. [jl4|, we demanded that the hardest two jets have energies ii^jet, , > 15 GeV, and that the 



total missing px exceeds 35 GeV; this removes most SM backgrounds. Fig. 4 shows the energy 
spectrum of the third (softest) jet after these cuts. 
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Figure 4: The spectrum of the softest jet produced in titl production at a future linear collider 
after the cuts described in the text, where the events have been forced into a 3-jet topology 
using the fcy clustering algorithm. The small shoulder at Ejet.^ = 15 GeV is caused by the cut 
on the energies of the two hardest jets. The dotted curve shows the result in the absence of 
hard gluon radiation, as predicted using fragmentation scheme II, while the other three curves 
show the total inclusive cross section. They differ in the fragmentation scheme and/or in the 
value of the fragmentation parameter e, as indicated. 

The dotted curve shows the contribution of events without hard gluon (defined through an 
IR cut-off i?™™ = 5 GeV), while the other three curves include contributions from both 2 — >■ 2 
and 2 — > 3 processes. We see that for Ejet^ > 15 GeV, the cross section is dominated by pertur- 
bative contributions. However, a significant dependence on the details of non-perturbative ti 
fragmentation remains even in this region. We saw in Fig. 1 that for a given value of e, scheme 
Ie predicts a softer stop "meson" spectrum, i.e. harder fragmentation jets, than scheme II 
does. It is therefore not surprising that scheme Ie predicts the third jet to be somewhat more 
energetic. However, some difference between the predictions from these two schemes remains 
even if the value of e is adjusted to reproduce the same average stop "meson" energy. After 
this adjustments the two schemes make very similar predictions for -Ejetg < 30 GeV, where 
most events still have an effective 2^2 kinematics. Nevertheless, our different methods 
to ensure overall energy-momentum conservation still give somewhat different predictions for 
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-^jctg > 30 GeV, where most events contain a quite energetic gluon. 

These differences are partly due to the fact that we force all events into a 3-jet structure, 
so that the gluon is often merged with (some of) the fragmentation products. This source of 
dependence on the fragmentation scheme can probably be reduced if a different jet definition 
is used. On the other hand, the two schemes also predict a somewhat different gluon energy 
spectrum even prior to any jet merging. We therefore expect that some dependence on the 
fragmentation scheme will remain for all jet definitions. Of course, for the parameters chosen 
in Fig. 4, the large number of tit* pair events (~ 65, 000 signal events after cuts for an inte- 
grated luminosity of 500 fb^^) should also allow to study various aspects of non-perturbative 
ti fragmentation in detail. Moreover, since the differences between the predictions of our 
fragmentation schemes are far smaller than the overall contribution of events containing hard 
gluons, even our current level of understanding of non-perturbative ti fragmentation would 
suffice to ascertain with a high degree of confidence that ti is indeed a color triplet. 



4) Summary and Conclusions 



We have examined the effects of both perturbative and non-perturbative strong interactions 
on titl production at current and future e~^e~ colliders. The latter manifests itself in ti 
fragmentation into a stop "meson", while the former gives rise to events containing a hard 
gluon in addition to the stop squarks. 

Non-perturbative ti fragmentation has been discussed in some detail in Sec. 2, with par- 
ticular emphasis on the conceptual ambiguities of the fragmentation of heavy partons. These 
ambiguities are largest for beam energies somewhat, but not far, above the pair production 
threshold, i.e. precisely in the region where squark production can best be studied at e~^e~ 
colliders. They include not only the definition of the fragmentation (scaling) variable, but also 
the algorithm used to ensure overall energy-momentum conservation. The scaling variable can 
be defined either through the energy or through the absolute value of the 3-momentum (or 
through some combination thereof). These two quantities coincide for massless partons, and 
also for massive partons in the infinite momentum frame where QCD factorization theorems 
are often formulated. However, they can differ greatly for massive partons in the laboratory 
frame, where event generators are needed to describe the data. We also introduced two dif- 
ferent schemes to guarantee overall energy and momentum conservation after fragmentation. 
The first scheme is based on the re-scaling of all 3-momenta, and is quite similar to algorithms 
employed in multi-purpose event generators like ISAJET. In the second scheme, which bears 
some resemblance to algorithms often used to model perturbative parton showering, one ex- 
plicitly generates a distribution in the virtuality of the produced squarks. We saw that these 
schemes can lead to quite different results when applied naively to the same fragmentation 
function. However, at least in the simple case of exclusive tit* production, i.e. in the absence 
of hard gluon radiation, most of these differences can be absorbed in a redefinition of the 
parameter(s) of the fragmentation function. 

Events containing hard gluons have been studied in Sec. 3. After giving explicit expressions 
for the corresponding fully differential cross section, we showed that such events should allow 
to extend the sensitivity of e~^e~~ collider experiments to scenarios with arbitrarily small mass 
splitting between ti and the lightest superparticle Xi? which we assume to be invisible. The 
gluon by itself can produce a sufficient amount of visible energy, and of missing transverse 
momentum and acoplanarity, to discriminate between signal and background even if ti decays 
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release almost no visible energy. We found that the cuts used to suppress Standard Model 
backgrounds also reduce the sensitivity of the predicted cross section to the details of ti 



fragmentation to a negligible level. Hard gluon radiation had previously been found [0 to 
increase the sensitivity of hadron colliders to tit^ pair production in scenarios with small 
stop-LSP mass splitting. However, that analysis treated gluon radiation only in leading 
logarithmic approximation (through the parton shower algorithm of ISA JET), rather than 
through exact matrix elements. Moreover, scenarios with mass splitting below 20 GeV still 
remained undetectable in hadronic collisions. 

We have also indicated how an analysis of 3-jet events from titl pair production at a future 
linear collider operating at high energy and high luminosity can be used to determine the color 
charge of ti. In this case some sensitivity to ti fragmentation remains after cuts. 

Our numerical results include both real and virtual ordinary QCD corrections. However, 
we did not include initial state radiation, nor did we include corrections from loops including 



gluinos, i.e. "SUSY QCD corrections" |10, 16 1. Both effects mostly change the overall nor- 



malization of the titl production cross section, leaving the shapes of the distributions studied 
here almost unchanged.^ The total cross section also depends on the ti — Ir mixing angle 9i, 
which we set to zero in our numerical examples. 

The numerical analyses presented here are by no means an exhaustive list of the effects 
of strong interactions on stop production at e^e~ colliders. For example, strong interactions 
will play a role in measurements of the stop and LSP masses from tit\ events. This is true in 
particular in scenarios with small stop-LSP mass splitting, where a large fraction of the total 
visible energy can come from non-perturbative ti fragmentation or hard gluon radiation, as 
illustrated in Fig. 2. Even if this mass splitting is quite large, the extracted values of these 
masses could be shifted by several GeV if the analysis is not designed such as to minimize 
these effects. 

While the other very short-lived squarks are not expected to undergo non-perturbative 
fragmentation, hard gluon radiation will occur with equal relative frequency for all squark 
flavors. The cross section given in Sec. 3 can easily be adapted to the production of other 
squarks by changing the relevant electroweak couplings. Many SUSY models predict first and 
second generation squarks to be closely degenerate in mass; moreover, SU{2) doublet squarks 
might undergo quite complicated cascade decays even if they are lighter than the gluino. One 
will then need a more sophisticated analysis to experimentally determine their color charges. 
The effect of hard gluon radiation on squark mass measurements in such a scenario []T1| also 



needs to be investigated. On the other hand, since ti might well be very long-lived compared 
to the time scale of strong interactions in spite of being very heavy in relation to the QCD 
scale A, it might even provide us with a unique tool to study non-perturbative aspects of 
the hadronization process. Analyses of squark pair production at e^e~ colliders can thus not 
only yield valuable information on these (as yet hypothetical) new particles, but might even 
improve our understanding of some basic aspects of strong interactions. 
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